arXiv:1501.05786v4 [astro-ph.CO] 8 Dec 2015 


Constraints on hybrid metric-Palatini models from background evolution 

Nelson A. LirngQ and Vanessa Smer-BarretcH] 

Institute for Astronomy, University of Edinburgh, 

Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK 

In this work, we introduce two models of the hybrid metric-Palatini theory of gravitation. We 
explore their background evolution, showing explicitly that one recovers standard General Relativity 
with an effective Cosmological Constant at late times. This happens because the Palatini Ricci 
scalar evolves towards and asymptotically settles at the minimum of its effective potential during 
cosmological evolution. We then use a combination of cosmic microwave background, supernovae 
and baryonic accoustic oscillations background data to constrain the models’ free parameters. For 
both models, we are able to constrain the maximum deviation from the gravitational constant G 
one can have at early times to be around 1%. 
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I. INTRODUCTION 

The ACDM model, which provides the simplest ex¬ 
planation for the Universe’s accelerated expansion in the 
form of a cosmological constant A (for a review on this 
topic, see m), is the main framework of modern cosmol¬ 
ogy. Despite being in agreement with supernovae ob¬ 
servations [ng, data from the cosmic microwave back¬ 
ground (CMB) [51 [7] including the recent Planck data 
[5] , and large-scale structure (LSS) data |5] , cosmologists 
still struggle to account for the difference between the 
theoretically-expected value for its energy density and 
the observed one. 

In light of these issues, a new physical model may be in 
order to account for that major component of our Uni¬ 
verse, usually labeled dark energy (DE). Theories that 
obey the standard gravitational field equations of Gen¬ 
eral Relativity (GR), such as quintessence, k-essence, and 
so on, propose scalar fields rolling in a potential (see [10] 
and references therein for a comprehensive review). How¬ 
ever, a much larger set of models arise from investigations 
where GR is assumed to fail, therefore proposing correc¬ 
tions to Einstein’s action. These are grouped under the 
name of Modihed Gravity Theories (MGT), such as the 
Brans-Dicke scalar-tensor theory El, higher dimensions 
like in braneworld models Enn], Galileon models [14], 
the Fab Four [TS] , f{R) theories [T5] . and many others. 
For an extensive review on MGT, see El- 

Of particular relevance to the background of the 
present work are f{R) theories. In these models, GR is 
modified by taking the Lagrangian density to be a non¬ 
linear function / of the Ricci scalar. The dynamics of 
this metric variational approach can be interpreted as 
due to a Brans-Dicke scalar field with parameter w = 0 
and a non-trivial potential V{(/)). However, when Solar 
System constraints are imposed, the interaction of this 
scalar held is restrained to a milimetric scale, nullifying 
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its impact on cosmology m- 

Another class of gravity theories, known as Palatini, 
emanate from considering the metric and the connection 
as independent variables m- However, these models 
have been observed to be unstable to matter loop correc¬ 
tions. If the matter loop effects are accounted for, they 
give rise to a Lagrangian density for the gravitational ac¬ 
tion which is a function of the Ricci scalars of both the 
connection and the metric, which in turn creates a set of 
theories named hybrid-metric Palatini [20l |2T|. The lat¬ 
ter are the focus of this paper, as we deem interesting to 
explore the late time acceleration effects of a set of mod¬ 
els who consider perturbative quantization methods |18j 
whilst also displaying connections with non-perturbative 
quantum geometries [20] . 

Like the pure metric and Palatini cases, the hybrid 
theory has a dynamically equivalent scalar-tensor repre¬ 
sentation [20l [21]. The authors of those papers have also 
shown that the scalar field need not be massive in or¬ 
der to pass the stringent Solar System constraints m, 
in contrast to the metric f{R) theories, while possibly 
modifying the cosmological [22] and Galactic [23] dynam¬ 
ics due to its light, long-range interacting nature. It has 
actually been shown that the additional scalar degree of 
freedom introduced by this theory does not need to be¬ 
come massive in order to evade solar-system tests. The 
viability of the theory on the smallest scales can be as¬ 
sured as long as the background value of the scalar field 
remains small [2Q1I2I]. Hence, this theory seems to not 
need an explicit screening mechanism, even though much 
work remains to be done on this topic. 

The criteria for obtaining cosmic acceleration have al¬ 
ready been discussed and introduced [22] ■ Alongside 
that, several cosmological solutions were derived, de¬ 
pending on the form of the effective scalar field potential, 
describing both accelerating and decelerating Universes. 
Lastly, the full set of linearly perturbed Einstein equa¬ 
tions was derived for this theory, and the evolution of the 
metric potentials was shown for a designer model repro¬ 
ducing the ACDM background evolution exactly [23]. 

In this work, we show the background history predicted 
by two /{TZ) models we introduce, which is achieved by 
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numerically evolving the Palatini Ricci scalar, TZ, that is 
then used to compute the remaining background quanti¬ 
ties. We then test the models against background data, 
using this to constrain their free parameters. Therefore, 
in Section [nj we briefly review the hybrid metric-Palatini 
theory formalism, focusing on the background equations 
we will be using throughout this work. In Sec. m we 
introduce and dissect the two models we have used in 
this study, showing the background evolution predicted 
by these against ACDM. In Section |IV| we present the 
background constraints on our models using a MCMC 
analysis, and finish in Sec. |V] with some concluding re¬ 
marks on this work. 


II. COSMOLOGY IN THE HYBRID 
METRIC-PALATINI GRAVITY 


The four-dimensional action describing the hybrid 
metric-Palatini gravity is given by 

^ ^ [R + /(7^)] + J d^x^C^, (1) 


where = 8nG and we set c = I. , V') is 

the standard minimally-coupled matter Lagrangian, R is 
the metric Einstein-Hilbert Ricci scalar and TZ = g^'^TZ^v 
is the Palatini curvature. The latter is defined in terms of 
the metric elements, g^'^, and a torsion-less independent 
connection, T, through 

^ = g^'' - f+ f. (2) 

Varying the action 0 with respect to the metric, one 
obtains the usual set of Einstein equations, given by 

G/ij/ + F{JZ)TZ^i, — -f{TZ)gfj.i^ = (3) 

where is Einstein’s tensor, and F{JZ) = df{JZ)/dTZ. 
is the matter field’s stress-energy tensor, defined as 


^ ^_2 6{^£^) 

Hgn ■ 

Equation 0 can be traced, yielding 

R + k^T = F{TZ)n - 2f{TZ), 


(4) 

(5) 


where T is the stress-energy tensor trace. The last equa¬ 
tion shows that the modifications to Einstein gravity are 
controlled by the failure of the standard GR trace equa¬ 
tion. 

On the other hand, varying the action Q with respect 
to the independent connection, r“j^, one gets the follow¬ 
ing equation, 

v„ {V^Fin)g^'') = 0, (6) 


which implies that is the Levi-Civita connection of 
a metric = F{TZ)g^i,. Therefore, is a conformal 


transformation of the original metric elements by the 
factor F{TZ). Hence, one can verify that the Palatini 
Ricci tensor, TZfj,^, is related to the metric one, R^y, by 
the following relation 

R-piv = Rp,u + 2 (TZ) ~ 

Adopting a flat Friedmann-Robertson-Walker (FRW) 
metric, = —dt'^ + a^{t)dx‘^, one can get the modi¬ 
fied Friedmann equations: 


3H^ 


and. 


1 

1 + /t?. 


k'^P - 3Hfn 


4/7^ 2 

( 8 ) 


2H = 


I 

1 + /t?. 


—{p +p) + FIJti — jn + 


ifn 

2-fn 


(9) 


where a dot stands for a differentiation with respect to 
time, t, and iJ = d/a is the Hubble parameter; p and p 
are the total energy density and pressure of the Universe. 
An equivalent equation for FI can be taken from the trace 
of Eq. yielding 

i/= 1 [F(7^)7^ - 2/(7^) - - 1277^] , (iQ) 

where we have used R = 6 (^H -|- 2H^^. 

A scalar-tensor description of this theory has been de¬ 
veloped, where the additional scalar degree of freedom is 
identified as </> = F(TZ) [501HJ. Eq. (II) in is an ef¬ 
fective Klein-Gordon equation for the scalar field. Using 
the FRW metric, it takes the form 




t 

2(/ 


0 

3 


2U - (1 -f </) 


dV 

d(j) 



T, (II) 


Equation ( [TT] ) can be re-expressed as a dynamical equa¬ 
tion for the Palatini Ricci scalar, TZ. Recalling that 
(j) = F{TZ) = Jti, then one can set F = TZfnn^ where 
fnn is the second derivative of fiJZ) with respect to TZ. 
A similar procedure can be done for higher order deriva¬ 
tives, allowing to rewrite Eq. 0 as 


7 ^ = -- 


I 


fnn 
fn 


fnnn 


_ fnn \ 
Vn) 

[TZ ifn - 1) - 2/(7^)] - K 


+ iHTZfnn + 

2 fn rp 

T 


( 12 ) 


where we have used Eq. (5) in [50], and funn is the 
third order derivative of f{R) with respect to R. From 
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this equation, we can define an effective potential where 
the effective Palatini Ricci scalar will roll, since 

^ ^ [7^ (/t^ - 1) - 2/(7^)], (13) 

meaning V(7Z) will be given by the indefinite integral 


P(7^) = 



[7^ (M - 1) - 2/(7^)] d7^. 


(14) 


It was shown that the hybrid metric-Palatini theory 
reduces to General Relativity with a possible cosmolog¬ 
ical constant in vacuum, since it shares the property of 
pure Palatini /(i?) theories in Minkowski fiat space-time 
|25j . Furthermore, the field’s equation of motion have 
been analyzed as a dynamical system: it was explicitly 
shown that as long as one provides a suitable V(^), or 
equivalently a function f(7^) such that the slope of the 
potential is downwards and its minimum happens for a 
small value of the scalar field (/>, one should always obtain 
a natural transition from standard cosmological evolution 
to accelerated expansion close to the present while also 
avoiding any conflict with solar-system constraints |22) . 


Hence, we use Eqs. (101 and ( |12[ ) to numerically evolve 
the background quantities predicted by a specific /{TV) 
model. To set the initial conditions, we have fixed a very 
small value for f-ji at a high redshift, z;, such that the de¬ 
viation from the Gravitational constant, G, is effectively 
small Eoiiiiiiiaiii] in the high curvature regime. Then, 
one can invert F to find TZ at that redshift. 


Lastly, we set T^i = 0 to minimize deviations from stan¬ 
dard General Relativity and use Eq. to solve for the 
Hubble parameter at Zi. The Ricci scalar, i?, can be com¬ 
puted using the GR relation, R = Q (h + . Even 


though T^i = 0 is a strong assumption, we have tested 
the models for a fairly broad range of initial velocities, 
within a slow-roll regime, and observed that their qual¬ 
itative behavior remained unaltered with the late-time 
evolution tending to an effective Gosmological Constant, 
and TZ asymptotically reaching the equilibrium position 
predicted by the effective potential we define in Eq. (14). 


III. MODELS OF HYBRID METRIC-PALATINI 
GRAVITY 


In this section, we introduce the models we defined 
for the hybrid metric-Palatini theory. While the gen¬ 
eral framework of this theory was derived in Refs. [2(H - 
I23j . they did not write down specific models and explore 
their consequences. These were inspired in theories of 
f{R) gravity, such as the Starobinsky [26] and the expo¬ 
nential m models, but are essentially phenomenological 
choices that are simply tested with the background evo¬ 
lution formalism displayed before. 


A. The exponential model 


The first model we introduce is defined by an exponen¬ 
tial function, given by 


/(7^) = A* (l + , 


(15) 


where A* and TZ^, are the model’s parameters, both of or¬ 
der Hq, where Hq is the present-day value of the Hubble 
parameter. We choose to define 7?.* as a positive con¬ 
stant, while A* should be a negative constant, since it 
should dominate over the other corrective terms intro¬ 
duced by our models. This is particularly relevant at 
late-times where, to recover standard GR with an effec¬ 
tive Cosmological Constant, one should have A* « —2A. 
Hence, A* < 0 allows the effective Cosmological Con¬ 
stant to have the correct sign, which becomes clear, for 
instance, looking at Eq. (§ in vacuum. 

The effective potential, V{TZ), associated to this model, 
is given by doing the indefinite integral defined by 
Eq. (14), which has the simple form 


V{TZ) = ^ 


2A*7^ -k ^TZ^ - 2A*iR,*e-^/^* 


We can see that the potential will be mostly dominated 
by its quadratic terms for order unit values of 7?.*. The 
minimum can then be estimated by 

TZnun « -2A*. (17) 

The effective cosmological constant value at which the 
Ricci scalar should sit in vacuum is given by the trace 
equation 


Aeff = i [-2 A*/k(-2A*) - 2/(A*)] « -A,. (18) 

Hence, we do expect to obtain a ACDM like evolution in 
the distant future, if our solution for TZ is to settle at the 
minimum of its potential. 

We choose a small value for f-j^ (less than unit), inde¬ 
pendently of the starting redshift Z;. The reason for this 
is, as mentioned in the end of Sec. jnj that f-ji sets the 
initial deviation from the actual gravitational constant 
in the high curvature regime. Therefore, higher values of 
f-ji will either suppress gravity considerably or invert its 
sign if it is smaller than —1. 

In Fig. we plot the background evolution predicted 
by this model for a set of parameters as a function of 
redshift, z, using the prescription described at the end of 
Sec. |Hj We use a Brent algorithm (see [55]) to find the 
correct A* value that recovers a flat cosmology. Hence, 
the only true free parameters will be 77.* and fn{zi)- 

As we can see in Fig. the evolution of the Palatini 
Ricci scalar starts at a position where the potential is 
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FIG. 1. The background evolution predicted by the expo¬ 
nential f{TZ) model compared to ACDM. We choose to plot 
TZ far into the future (2 —> —1), to explicitly show that our so¬ 
lution asymptotically tends to the minimum of the potential, 
V{TZ), which we plot as well. We also plot the deceleration 
parameter q. The present-day matter energy density, flm, is 
set to 0.30. 7?.* is in units of Hq. 


tilted and one would expect for it to roll down towards 
the minimum. However, at early times, the evolution 
of TZ is dominated by the stress-energy tensor trace or, 
equivalently, by the matter energy density. Since the ra¬ 
tio fn/fnn is negative throughout the whole evolution, 
the matter density contribution pushes TZ upwards the 
effective potential, while its slope and the Hubble fric¬ 
tion term exert the opposite effect. As matter ceases 
to dominate close to the present, TZ inverts its motion 
and starts evolving towards the minimum, where it will 
asymptotically settle in the distant future. 

Lastly, in Fig.j^ we analyze the deceleration parameter 
q. We observe that our model predicts a Universe that 
will be expanding in an accelerated manner today, tran¬ 
sitioning from a matter dominated decelerating phase at 
a redshift of around z « 1, as we have g < 0 at z = 0. 
This is a general result of this model, as can be inferred 
from the deceleration parameter equation 


q = 


H 

IP 


- 1 , 


(19) 


where H is given by Eq. 10 Switching from physical time 


t to In a we will have, at z = 0, 


q = - 


ml 


[F{TZ)TZ - 2/(7^) -b - l2Hl] - 1. 

( 20 ) 

As we observe in Fig{T] our model evolves towards small 
values of F{TZ) today, as the exponential is suppressed 
by the larger values of TZ. Therefore, in Eq. B we can 
neglect the F{TZ)TZ term. Then, from Eq. 15 we note 
that f{TZ) will be dominated by A*. Since this parameter 
is determined from imposing a flat cosmology, A,^ « ~2A, 
where A « 2,.\Hl is the actual cosmological constant. In 
light of these arguments, Eq. [^becomes 


We verify in Fig. that our prediction matches remark¬ 
ably well the numerical values of obtained for q today for 
different parameters of the exponential model. Hence, 
we conclude that this model should always predict an 
accelerated expansion today. 


B. The quadratic model 


The second model we introduce is the quadratic model, 
which we define by the function 


/(7^) = A* (1 + 7^V7^*) 


( 22 ) 


where A* and 7?.* are a negative and positive constant of 
order TLg, just like in the previous model. 


Computing the indefinite integral defined in Eq. (14) 


one can find the effective potential associated to this 
model: 

V(TZ) = -^-k.TZ^-^-TZ\ 


(23) 
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which clearly has a global minimum at 72, = 0 and a 
maximum at 72. = 2A*. Therefore, we expect the solution 
for 72, to asymptotically settle at the minimum in vacuum, 
leading to an effective cosmological constant value of 

Aeff = i[-2A,/7^(0)-2/(0)] =-A*. (24) 

In Fig. we can see the evolution predicted by the 
quadratic model for H, compared to ACDM’s, for differ¬ 
ent values of fn{zi) and 72.*. We have chosen a negative 
fn{zi) as this sets the initial TZ to be positive. 

In Fig. we observe that 72 rolls down the effective 
potential towards the minimum from the beginning of 
the evolution, asymptotically settling there. When com¬ 
pared to the behavior seen for the exponential model 
in Fig. the different dynamics are linked to the ra¬ 
tio fn/fun multiplying the stress-energy tensor in the 
dynamical equation for 72.. For this particular setting of 
the quadratic model, this ratio is positive at early times, 
meaning that 72 will be pulled down by the matter en¬ 
ergy density term towards the minimum of the potential, 
asymptotically settling there when the matter contribu¬ 
tion becomes negligible. Had we chosen 72.(zi) to be posi¬ 
tive instead, the behavior wouldn’t qualitatively change, 
with the matter term dictating 72. to evolve towards the 
minimum in a symmetric manner, starting from 72 < 0. 

Also, note that by setting the initial conditions for 
fnizi) such that deviations from standard GR are kept 
small, 72. starts its evolution already close to the mini¬ 
mum, guaranteeing that it settles onto the minimum of 
the effective potential V{TZ). Had we started the evolu¬ 
tion such that the initial value of 72. would be beyond the 
maximum observed for the effective potential in Fig.[^ 72, 
would roll indefinitely, jeopardizing the late-time achieve¬ 
ment of an effective cosmological constant in this model. 

Lastly, we have the evolution of the deceleration pa¬ 
rameter q in Fig. [^predicted by this model. We observe 
that, close to the present, /(72) will be dominated by A* 
as can be seen in Eq. |22[ since 72 is tending to the min¬ 
imum of the potential V (72) located at zero. Hence, the 
analysis of the parameter q in Sec. |HI A| applies in this 
case as well. Therefore, this model should always predict 
an accelerating Universe today. 


IV. BACKGROUND CONSTRAINTS 

A. Observables 

To constrain our models parameters, we perform a 
Metropolis-Hastings analysis using several background- 
related observables: the luminosity distance from Su¬ 
pernova type la (SNIa); the baryon acoustic oscillations 
(BAG); the shift parameter, acoustic scale and redshift 
of decoupling from the cosmic microwave background 
(CMB). For a detailed explanation on the statistics treat¬ 
ment of SNIa, BAG and CMB see ^5] . 





FIG. 2. The same as Fig. for the quadratic f(TZ) model. 
We have again fixed Um = 0.30. 


1. SNIa luminosity-redshift relation 

For the supernovae analysis, we have used the Union 
2.1 SNIa catalog from the “Supernova Cosmology 
Project” (SCP) [SD]- This data set is a compilation of 
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TABLE I. Inverse covariance matrix for the distance informa¬ 
tion obtained from Planck in the ACDM framework. 


Planck 



Best fit 

Ia{z*) 

R(z*) 

21 * 

Ia{z*) 

301.77 

44.077 

-383.927 

-1.941 

R{z,) 

1.7477 

-383.927 48976.330 

-630.791 

2 :* 

1090.25 

-1.941 

-630.791 

12.592 


580 type la Supernovae located over the redshift interval 
0.623 < z < 1.415. In our Metropolis-Hastings analysis, 
we have marginalized over the nuisance parameter Hq 
using the procedure described in the appendix [5T| . 


8. Baryon Acoustic Oscillations peak 

Following the Planck analysis [32], we have used 
BAO observations from the 6dF Galaxy Redshift Sur¬ 
vey (6dFGRS) at low redshift, = 0.106) = 0.336 ± 

0.015 [32, from the Luminous Red Galaxies (LRG) sam¬ 
ple of the Sloand Digital Sky Survey (SDSS) 7-year 
data release at the median redshift, d°^^{z = 0.35) = 
0.1126 ± 0.0022 [33]) and the BAO signal from BOSS 
CMASS DR9, d°^%z = 0.57) = 0.0732 ± 0.0012 [35|. 


3. Cosmic Microwave Background 

In this work, we use the Planck distance information 
for the AGDM model to constrain our models. In table 
m it is possible to see the inverse covariance matrix and 
best-fit values obtained from Planck [33] - 


B. Metropolis - Hastings Algorithm 


we have chosen to impose over them. Our two /{TZ) 
models have a set of three parameters they share. These 
are the present-day relative energy densities of matter 
and baryons, and fib, and the present-day value of 
the Hubble parameter, Hq. The ranges we have chosen 
for them are, respectively, [0.01,0.99], [0.001,0.080] and 
[40.0,100.0]. 

Also for both models, A* is determined by the back¬ 
ground evolution, assuming we recover a flat cosmology 
today, at a = 1. As described in Sec. [Ill A[ this is 
achieved using a Brent algorithm (see m) which en¬ 
sures that the present-day value of the numerical Hub¬ 
ble parameter we obtain coincides with Hq. We start 
our background evolution at a redshift of « 1 x 10®, 
and neglect any relativistic effects. Hence, we take the 
equation of state for matter to be zero throughout, i.e. 
Pm oc a“®. 

For the exponential model, we study two situations. In 
the first case, we choose to keep fn at Zi fixed to a very 
small value, 10“"^, and let 77.* vary between [1.0,10.0]. 
We do so to test the data against a definitive modifica¬ 
tion of gravity where the AGDM limit is not explicitly 
attainable. 

In the second case, we let /^(zi) vary between [1 x 
10“®,0.1] and 77* between [0.01,15.0]. We can not con¬ 
sider fnizi) = 0 due to the way we set the initial condi¬ 
tions in this model, as this would lead to a logarithmic 
divergence. Therefore, we can not set the deviations from 
standard GR plus ACDM at early times to zero, but only 
asymptotically minimize them by taking /^(zi) —>■ 0^. 
Nevertheless, we can test how large a deviation should 
be possible by considering an upper limit for of 

order 1 x 10“^. 

Lastly, for the quadratic model, we fix 77* to a chosen 
value and let f-ji at z^ vary between [—0.1,0.lj. We chose 
to fix 77* since the effective potential on which 77 evolves 
is independent of 77* in the quadratic model, as seen in 
Eq. ([2^. 


Having described the procedure to calculate the 
of each observable we consider in the previous section, 
we continue towards the calculation of the confidence 
contours by means of a Metropolis-Hastings algorithm, 
which is a Markov Chain Monte Carlo method based on 
a stochastic sampling technique m- One of the main 
advantages of the Metropolis-Hastings algorithm is its 
treatment of marginalized variables. When considering a 
subset of the parameters that form a chain, the marginal¬ 
ization over the remainder of the parameters occurs im¬ 
mediately, therefore making the treatment of the chains 
a simple process 


C. Priors 

In this section we discuss the range of variation of the 
parameters we use in our models, which are the flat priors 


D. Results 


In this section, we present the marginalized 2-d con¬ 
tours for the posterior probabilities distributions of our 
parameters which were calculated using a Metropolis- 
Hastings algorithm, as described in Sec. jlVBj The plots 
shown here exhibit the combined constraints of the three 


background surveys we described in Sec. jlV A| These 
plots were done using the plotting functions available in 
CosmoMC |35] - 

In Fig. [^ we have a triangular plot with the 2-d con¬ 
tours between TLq, Dm and 77*, as well with the indi¬ 
vidualized posterior distributions of each parameter for 
the exponential model considering fnizi) fixed to 10“"^. 
We obtain a slightly higher value of Hq than the recent 
Planck result [32], but we do not perform such a com¬ 
plete analysis, limiting ourselves to background observ¬ 
ables. We also observe a smaller Dm value than in 1321. 
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^0 fni^i) 


FIG. 3. The 2-d contours of the combined constraints from 
the background surveys we consider for the exponential f{TZ) 
model with fixed to lO”"*. We also plot the individ¬ 

ual marginalized posterior probability distributions of each 
parameter. 



FIG. 4. The 2-d contours of the combined constraints from 
the background surveys we consider for the exponential f{TZ) 
model by considering a prior range for fiz{zi) between [1 x 
10“®, 0.1]. We also plot the individual marginalized posterior 
probability distributions of each parameter. 

which results from the combination of the different sur¬ 
veys we have considered, as the Union2.1 and BAO sur- 


FIG. 5. Similar as Figs. and for the quadratic f{TZ) 
model. 


veys do tend to prefer a slightly smaller flm value than 
Planck [30l [32] . The l-cr limits on these parameters are 
Ho = 68.9 ± 0.7 and Qm = 0.303 ± 0.009. 

For the 77.* parameter we cannot clearly state the con¬ 
fidence limits, as these are completely prior determined. 
Interestingly, we do observe a preference towards smaller 
values of 77.*, possibly extending all the way to 0 had 
we considered that limit. We set our initial conditions 
by imposing an initial f-jz value at the starting redshift Zi 
that we invert to obtain the corresponding 77i. Hence, for 
this model, 77i = —77.* In [—/7^(zi)77.*/A*]. Therefore, we 
cannot have a pure GR plus ACDM for the exponential 
model because neither /^(zi) or 77.* can be set exactly to 
zero. However, the closest this model can get to ACDM, 
for a fixed /^(zi), is when 77* —>■ 0+: in this limit, we 
observe that TZ\ tends to decreasingly smaller values as 
77* —?> 0+, while keeping the 'R-j'R.-i, ratio considerably 
large such that /(77) —?► A*. Hence, we recover an almost 
ACDM like evolution, which can be understood looking 
at the trace equation, Eq. ([^, which tends increasingly 
closer to the GR plus ACDM limit of i? -I- k^T = 4A. 

Fig. [^presents the combined constraints on the expo¬ 
nential model considering the flat prior on f'jzizi) between 
[1 X 10“®,0.1]. We can observe a clear upper limit on 
fizizi) of around 10“^, limiting, therefore, the maximum 
deviation from the actual gravitational constant G we 
can have at early times. However, now 77* appears even 
more unconstrained by the data, as larger values of 77* 
are also allowed since these can be compensated by the 
fnizi) —>■ 0"^ values: this limit pushes the model closer to 
the standard GR plus ACDM limit even for large values 
of 77.*. We reinforce that we can not set the deviation 































from standard GR exactly to zero because that would 
lead to a logarithmic divergence in the initial conditions 
of this model. However, the lower limit we have chosen 
for the prior is much smaller than the current pre¬ 

cision with which one can measure early deviations from 
the gravitational constant [39l EQI • 

Lastly, in Fig. we have the triangular plot with the 
confidence contours and 1-d marginalized posterior prob¬ 
ability distributions for the f{TZ) quadratic model. We 
obtain very similar results for the standard cosmological 
parameters Hq and flm as those observed in the expo¬ 
nential model in both situations, with the l-cr limits on 
them being Hq = 68.9 ± 0.8 and = 0.303 ± 0.010. 

For the quadratic model, we chose to keep 7^* fixed 
because it would not have a significant impact on the 
background evolution predicted by this model since it 
does not alter the shape of the potential on which TZ 
evolves. Hence, the third parameter we show constraints 
for is the initial value we impose for /t^ at the starting 
redshift, Zi. This, as detailed before, sets the maximum 
deviation from the standard gravitational constant G one 
can have at early times in this model, since f-ji evolves 
asymptotically to 0 from its starting value. 

It is clear that as fnizi) 0 we get closer to a ACDM 
like evolution. Given it would be numerically hard to 
evolve the model if fnizi) = 0, we made the approxima¬ 
tion that, if the Metropolis-Hastings algorithm encoun¬ 
tered such a value, we would have an exact AGDM back¬ 
ground evolution. The results we have obtained show a 
preference for a standard GR plus AGDM scenario, as 
can be seen in the 1-d posterior probability distribution 
for fTi{zi) in Fig. The corresponding l-cr confidence 
limits are fnizi) = —0.001 ± 0.006. We also observe a 
symmetry on the posterior distribution of this parame¬ 
ter, which could be expected given that the evolution of 
TZ is symmetric under the change of sign of f-ji for, as 
discussed in Sec. ElBl 


V. CONCLUSION 

In this work, we explored a way to obtain the back¬ 
ground evolution for two different models of the novel 
hybrid metric-Palatini theory of gravity. We re-wrote 
the dynamical equation for the additional degree of free¬ 
dom introduced by this theory as a dynamical equation 
for the actual Palatini Ricci scalar TZ. We define the ini¬ 
tial conditions by imposing the deviation one has from 
standard GR at early times. Hence, we set a small value 
for fnizi) and invert the latter to obtain TZ{zi), while 
keeping 'JZ{zi) = 0. 

We define an effective potential V(TZ) where the Pala¬ 
tini Ricci scalar evolves and, if a minimum exists, TZ 
should asymptotically settle there in vacuum, so that 
one recovers standard GR plus an effective cosmological 
constant at late times. V{TZ) could potentially have a 
complicated form. However, for the models we introduce 
here, that is not the case. 


We present the exponential and quadratic f{TZ) models 
and show that the background evolution predicted by 
them does not deviate much from ACDM. This could be 
different, of course, had we decided to set the deviation 
from the gravitational constant G in the high-redshift 
regime to be large. Also, we explicitly show the effective 
potential V{TZ) for both models and TZ asymptotically 
tending to its minimum at late times. This is less obvious 
in the exponential model as the matter term in Eq. (12) 
initially drives the Palatini Ricci scalar up the potential, 
only for it to later slowly fall down towards the minimum 
due to the potential slope. 

We also study and present the evolution of the decel¬ 
eration parameter q for our models. We verify, numer¬ 
ically, that they predict an accelerated expansion today 
for different values of their parameters. Furthermore, we 
perform a simple analytical analysis of q and conclude 
that our models yield a present-day value for it around 
— 1/2, meaning they are suitable candidates for produc¬ 
ing cosmological acceleration today. 

We then use background CMB, BAG and Supernovae 
data to constrain the models. Keeping fuizi) fixed to 
10“'^ for the exponential model, we cannot state an ac¬ 
tual constraint on the 7?.* parameter. However, we do 
note how the data, as expected, seems to tend towards 
the AGDM limit. We believe that, had we not chosen to 
restrict the prior range in order to have a definitive mod¬ 
ification of gravity without the ACDM limit, we would 
see the lower range of our confidence contours in Fig. 
tending to 0+ in TZ^,. 

Still for the exponential model, when we impose a flat 
prior on between [1 x 10“®, 0.1] we observe a clear 

upper limit of order 10“^. This value marks the max¬ 
imum deviation one can have at early times from the 
actual gravitational constant G. Also, the data exhibits 
a marked tendency towards the lower limit of as 

fizizi) —>■ O'*" minimizes the early time deviations from 
the gravitational constant G, allowing the model to get 
asymptotically closer to the standard GR-I-ACDM limit. 

For the quadratic model, we constrained the initial 
value of f-fi while keeping 7?.* fixed which means that, 
effectively, we are again constraining the maximum de¬ 
viation one has from the gravitational constant at early 
times. As expected, the results indicate a preference to¬ 
wards no deviation at all, as in the standard GR plus A 
limit. We obtain fnizi) = —0.001 ± 0.006 as the confi¬ 
dence limits for our parameter. 

Hence, we see in this work that fnizi) could play an 
important role in the constraining of this theory, as it sets 
the deviation one observes from standard GR at early 
times. The background data we have used permits a de¬ 
viation around 1% from the actual gravitational constant 
G at early times, which is within constraints on Ges/G 
coming from BBN and GMB data [39l |40] . 

This, combined with the fact that the Newtonian po¬ 
tentials also exhibit a departing behavior from AGDM 
at early times [Mj, suggests that it would be very inter¬ 
esting to constrain these models using the latest Planck 
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